2024_anticoagulation
Aaron Mittel
2025-11-28
Aims and Hypotheses
- Aim 1: Determine the prevalence and differentiating features of patients with versus without aPTT-anti-Xa discordance among patients receiving ECMO at CUIMC in 2024.
Hypothesis 1: Discordance is a relatively common, yet under appreciated phenomenon, and varies by geographic location of care and patient attributes including demographics, severity of illness, and underlying disease process.
- Aim 2: Determine the implications of discordance between aPTT and anti-Xa on hemorrhagic and thrombotic complications in a cohort of patients who received ECMO at CUIMC in 2024.
Hypothesis 2: Discordance between aPTT and anti-Xa indicates the presence of coagulation system pathophysiology which is a risk factor for development of hemorrhage and thrombosis.
Classic systemic anticoagulation goals, anti-Xa vs aPTT, at NYP:
- 0.3 - 0.7 vs 51 - 83
Typical anti-Xa goal vs aPTT goals among heparinized ECMO patients:
0.1 - 0.3 vs 35 - 50 (LOW GOAL)
0.2 - 0.5 vs 45 - 65 (CLASSIC GOAL)
Methods:
Only first ECMO run used.
ECMO duration: Dr. Connie Nguyen manually reviewed heparin-exposed patients’ EMRs in June 2025 to identify true ECMO time off, as initial data review identified some labs that were drawn while on heparin but after the patient was off ECMO. The variable “time_off” in the original data was obtained from ELSO sources. Per Dr. Nguyen’s review, if her evaluation of ECMO stop time was within 4 hours of ELSO-derived time_off, then the actual true time off was deemed as that known by ELSO. If there was a greater than 4 hour discrepancy, a new time_off_true value was recorded.
Lab cutoffs:
if aPTT <20.0 - made 0
if aPTT >180 - made 180
if anti-Xa <0.04 or <0.1 - made 0
if anti-Xa >2.00 - made 2.0
if data is duplicated at a specific timepoint (e.g. 2 lab samples, exact same result time but different values), then taking the mean of the two exact same-timed values
Hyperbilirubinemia can artifactually elevate anti-Xa (colorimetric assay). Jaundice threshold is 2 - 3 mg/dL total bilirubin. Upper limit of normal total bilirubin at NYP-CUIMC is 1.2 mg/dL.
Bilirubin is measured in separate tubes compared to aPTT, anti-Xa,
INR, and fibrinogen. As hyperbilirubinemia can impact assessment of
anti-Xa but also may not change rapidly from one moment to another, we
accounted for total bilirubin by computing the
nearest_bilirubin value within 24 hours of any aPTT. (aPTT
chosen as it was the most frequently measured coagulation lab).
COMPLICATIONS: all complications reviewed were documented during ECMO duration, with exception of clinically significant bleeding and thrombotic complications which were manually assessed by Connie and Phil as follows:
DVT/PE_1_time_and_date
DVT/PE_2_time_and_date
DVT/PE_3_time_and_date
DVT/PE_4_time_and_date
Ischemic_stroke_1_time_and_date
Ischemic_stroke_2_time_and_date
Ischemic_stroke_3_time_and_date
Acute_limb_ischemia_1_time_and_date
Acute_limb_ischemia_2_time_and_date
Acute_limb_ischemia_3_time_and_date
GI bleeding - present on CTA and confirmed via endoscopy
ICH - based on imaging findings
Basic Lab Counts and Scatterplot visualization among ECMO paients who received heparin exposed
| n |
|---|
| 165 |
| department_name | n |
|---|---|
| MIL 5 CTICU | 76 |
| MIL 5 CCU | 41 |
| MIL 4 MICU A | 31 |
| HRT CARDIAC CARE | 25 |
| MIL OPERATING ROOM | 5 |
| MIL 4 SICU | 1 |
| MIL CARDIAC CATH | 1 |
| GBG 4 W CT ICU | 1 |
| n |
|---|
| 113 |
| department_name | n |
|---|---|
| MIL 5 CTICU | 56 |
| MIL 5 CCU | 31 |
| HRT CARDIAC CARE | 22 |
| MIL 4 MICU A | 6 |
| MIL 4 SICU | 1 |
| MIL OPERATING ROOM | 1 |
| MIL CARDIAC CATH | 1 |
| GBG 4 W CT ICU | 1 |
## # A tibble: 1 × 1
## correlation
## <dbl>
## 1 0.740Polynomial Fitting
This fits a second degree polynomial on aPTT and accounts for repeated measures per patient (i.e. multiple labs per individual patient) by including a random intercept for each patient. This is important as the repeated measures per patient are not independent - the labs are dependent on the baseline values (and the baseline anti-Xa between patients who have the same aPTT may not be the same, or vice versa).
95% prediction interval chosen to balance sensitivity and specificity.
|
|
## Adding missing grouping variables: `mrn`, `lab_result_time`| mrn | lab_result_time | activated_partial_thromboplastin_time | heparin_assay_quantitative |
|---|---|---|---|
| 1001491594 | 2024-08-24 22:45:00 | 96.9 | 0.36 |
| 1002467250 | 2024-12-12 18:30:00 | 63.5 | 0.26 |
| 1006905710 | 2024-10-10 18:02:00 | 71.3 | 0.07 |
| 1006905710 | 2024-10-11 06:23:00 | 56.2 | 0.16 |
| 1007129753 | 2024-04-12 06:13:00 | 43.5 | 0.10 |
| 1007129753 | 2024-04-12 18:41:00 | 42.9 | 0.10 |
| 1007129753 | 2024-04-13 06:08:00 | 43.8 | 0.10 |
| 1007129753 | 2024-04-13 18:00:00 | 45.2 | 0.10 |
| 1009477257 | 2024-12-02 15:04:00 | 41.9 | 0.09 |
| 1010147999 | 2024-02-20 18:35:00 | 60.3 | 0.10 |
| 1010147999 | 2024-02-21 12:35:00 | 52.6 | 0.10 |
| 1011043983 | 2024-06-17 17:05:00 | 75.4 | 0.33 |
| 1011043983 | 2024-06-17 22:39:00 | 93.2 | 0.37 |
| 1011043983 | 2024-06-18 04:34:00 | 93.8 | 0.36 |
| 1011043983 | 2024-06-18 09:40:00 | 106.5 | 0.40 |
| 1011043983 | 2024-06-18 16:00:00 | 103.4 | 0.43 |
| 1011043983 | 2024-06-18 22:46:00 | 89.5 | 0.37 |
| 1011043983 | 2024-06-19 05:51:00 | 106.0 | 0.43 |
| 1011043983 | 2024-06-19 11:46:00 | 122.7 | 0.46 |
| 1011043983 | 2024-06-20 11:59:00 | 104.6 | 0.48 |
| 1011077110 | 2024-05-20 06:06:00 | 58.1 | 0.14 |
| 1011077110 | 2024-05-21 06:42:00 | 70.1 | 0.23 |
| 1100547209 | 2024-05-02 14:06:00 | 119.4 | 0.60 |
| 1102078503 | 2025-01-07 10:13:00 | 36.9 | 0.04 |
| 1201199342 | 2024-09-26 12:28:00 | 76.7 | 0.34 |
| 1400032239 | 2024-03-18 23:46:00 | 45.6 | 0.10 |
| 1400969240 | 2024-10-13 23:20:00 | 39.7 | 0.05 |
| 1400969240 | 2024-10-14 18:13:00 | 38.5 | 0.06 |
| 1401183967 | 2024-03-30 18:34:00 | 102.0 | 0.30 |
| 1401327774 | 2024-11-03 03:56:00 | 122.1 | 0.38 |
| 1401327774 | 2024-11-03 11:13:00 | 94.3 | 0.38 |
| 1401327774 | 2024-11-03 19:02:00 | 63.2 | 0.22 |
| 1401468984 | 2024-02-24 07:04:00 | 42.8 | 0.10 |
| 1401468984 | 2024-02-27 18:47:00 | 55.3 | 0.20 |
| 1401639110 | 2024-08-20 06:59:00 | 64.4 | 0.16 |
| 1401766288 | 2024-07-31 06:09:00 | 60.3 | 0.21 |
## # A tibble: 1 × 2
## mean_aptt mean_antiXa
## <dbl> <dbl>
## 1 72.3 0.241| concordance_classic_simple | n |
|---|---|
| concordant | 13 |
| discordant aPTT above antiXa | 22 |
| high aPTT low antiXa | 1 |
## Adding missing grouping variables: `mrn`, `lab_result_time`| concordance_low_simple | n |
|---|---|
| concordant | 20 |
| discordant aPTT above antiXa | 15 |
| high aPTT low antiXa | 1 |
| mrn | lab_result_time | activated_partial_thromboplastin_time | heparin_assay_quantitative |
|---|---|---|---|
| 1001491594 | 2024-08-20 04:09:00 | 41.8 | 0.47 |
| 1001491594 | 2024-08-20 10:16:00 | 55.3 | 0.59 |
| 1001491594 | 2024-08-20 21:41:00 | 61.2 | 0.67 |
| 1001491594 | 2024-08-21 05:17:00 | 66.8 | 0.74 |
| 1001491594 | 2024-08-22 09:33:00 | 46.4 | 0.43 |
| 1003642179 | 2024-01-12 05:15:00 | 63.9 | 0.56 |
| 1009544965 | 2024-01-31 14:58:00 | 98.6 | 0.90 |
| 1009544965 | 2024-02-01 01:05:00 | 92.9 | 0.80 |
| 1009544965 | 2024-02-01 17:11:00 | 68.1 | 0.60 |
| 1009544965 | 2024-02-02 00:55:00 | 57.2 | 0.50 |
| 1010696012 | 2024-04-01 03:44:00 | 44.3 | 0.40 |
| 1010696012 | 2024-04-01 09:37:00 | 46.2 | 0.40 |
| 1010696012 | 2024-04-03 06:25:00 | 35.1 | 0.30 |
| 1101397979 | 2024-02-23 22:03:00 | 170.7 | 1.70 |
| 1101397979 | 2024-02-28 12:44:00 | 55.8 | 0.50 |
| 1102216636 | 2024-08-29 04:04:00 | 67.8 | 0.62 |
| 1201555171 | 2024-10-05 01:13:00 | 44.3 | 0.40 |
| 1201555171 | 2024-10-05 12:40:00 | 56.2 | 0.63 |
| 1201555171 | 2024-10-05 19:00:00 | 53.4 | 0.57 |
| 1201555171 | 2024-10-06 01:07:00 | 44.7 | 0.45 |
| 1201555171 | 2024-10-06 06:29:00 | 38.2 | 0.34 |
| 1401554006 | 2024-01-15 03:42:00 | 39.2 | 0.33 |
| 1401554006 | 2024-01-15 18:16:00 | 41.4 | 0.39 |
| 1401554006 | 2024-01-16 05:37:00 | 43.8 | 0.39 |
| 1401591974 | 2024-02-08 03:45:00 | 70.3 | 0.70 |
| 1401591974 | 2024-02-10 17:55:00 | 44.5 | 0.40 |
| 1401591974 | 2024-02-12 05:09:00 | 53.4 | 0.50 |
| 1401611898 | 2024-04-10 12:45:00 | 48.0 | 0.40 |
| 1401611898 | 2024-04-11 07:22:00 | 35.4 | 0.30 |
| 1401611898 | 2024-04-11 20:24:00 | 46.9 | 0.40 |
| 1401611898 | 2024-04-12 05:14:00 | 41.9 | 0.40 |
| 1401611898 | 2024-04-12 12:44:00 | 43.7 | 0.40 |
| 1401682296 | 2024-04-17 04:01:00 | 36.3 | 0.30 |
| 1401683148 | 2024-05-24 05:06:00 | 51.2 | 0.40 |
| 1401689750 | 2024-04-22 05:22:00 | 60.2 | 0.60 |
| 1401689750 | 2024-04-23 05:38:00 | 45.7 | 0.50 |
| 1401689750 | 2024-04-24 04:17:00 | 39.8 | 0.40 |
| 1401689750 | 2024-04-25 04:24:00 | 39.4 | 0.40 |
| 1401696145 | 2024-05-07 19:31:00 | 47.0 | 0.40 |
| 1401697961 | 2024-04-29 19:21:00 | 33.8 | 0.30 |
| 1401707171 | 2024-05-10 19:23:00 | 39.4 | 0.50 |
| 1401707171 | 2024-05-11 04:21:00 | 34.0 | 0.30 |
| 1401707171 | 2024-05-11 11:42:00 | 30.7 | 0.30 |
| 1401707171 | 2024-05-12 07:19:00 | 34.4 | 0.40 |
| 1401707171 | 2024-05-12 12:43:00 | 34.9 | 0.30 |
| 1401707171 | 2024-05-12 22:13:00 | 35.1 | 0.40 |
| 1401707171 | 2024-05-13 03:56:00 | 35.5 | 0.40 |
| 1401707171 | 2024-05-13 21:40:00 | 28.6 | 0.30 |
| 1401796928 | 2024-07-25 18:01:00 | 43.9 | 0.51 |
| 1401848048 | 2024-09-10 01:40:00 | 48.4 | 0.43 |
| 1401848048 | 2024-09-10 06:43:00 | 53.5 | 0.50 |
| 1401848048 | 2024-09-10 12:27:00 | 53.2 | 0.54 |
| 1401991907 | 2024-12-09 06:48:00 | 39.2 | 0.34 |
| 1401991907 | 2024-12-09 19:11:00 | 40.2 | 0.36 |
| 1401991907 | 2024-12-10 00:53:00 | 40.0 | 0.37 |
| 1401991907 | 2024-12-10 06:26:00 | 36.4 | 0.35 |
| 1401991907 | 2024-12-10 13:57:00 | 39.8 | 0.34 |
| 1401991907 | 2024-12-10 18:34:00 | 38.6 | 0.38 |
## # A tibble: 1 × 2
## mean_aptt mean_antiXa
## <dbl> <dbl>
## 1 49.6 0.474## Adding missing grouping variables: `mrn`, `lab_result_time`| concordance_classic_simple | n |
|---|---|
| concordant | 19 |
| discordant aPTT below antiXa | 38 |
| low aPTT high antiXa | 1 |
| concordance_low_simple | n |
|---|---|
| concordant | 22 |
| discordant aPTT below antiXa | 35 |
| low aPTT high antiXa | 1 |
| mrn | n |
|---|---|
| 1001491594 | 2 |
## # A tibble: 7 × 3
## # Groups: lab_result_time, unusual_type [7]
## lab_result_time unusual_type n
## <dttm> <fct> <int>
## 1 2024-08-20 04:09:00 aPTT < antiXa 1
## 2 2024-08-20 10:16:00 aPTT < antiXa 1
## 3 2024-08-20 21:41:00 aPTT < antiXa 1
## 4 2024-08-21 05:17:00 aPTT < antiXa 1
## 5 2024-08-22 09:33:00 aPTT < antiXa 1
## 6 2024-08-23 05:28:00 concordant 1
## 7 2024-08-24 22:45:00 aPTT > antiXa 1Study ID 5 (MRN 1001491594) had both discordance types (aPTT > antiXa and aPTT < antiXa).
Perhaps it would be worth retaining this patient in the model, as discordance could contribute to hemorrhage and/or thrombosis in separate instances.
| unusual_type | n |
|---|---|
| concordant | 477 |
| aPTT > antiXa | 36 |
| aPTT < antiXa | 58 |
## # A tibble: 6 × 4
## conc high low n
## <lgl> <lgl> <lgl> <int>
## 1 FALSE FALSE TRUE 5
## 2 FALSE TRUE FALSE 3
## 3 TRUE FALSE FALSE 77
## 4 TRUE FALSE TRUE 13
## 5 TRUE TRUE FALSE 14
## 6 TRUE TRUE TRUE 1## # A tibble: 3 × 2
## first_discord_type n
## <chr> <int>
## 1 aPTT < antiXa 19
## 2 aPTT > antiXa 17
## 3 concordant 77Interpreting:
77 patients: always concordant
17 patients: ever high only (aPTT > antiXa)
19 patients: ever low only (aPTT < antiXa)
1 patient: both high AND low
Patient-Level Demographics, based on development of any discordance
Note that one patient with both discordance groups is included in this analysis. This patient’s first discordance type was “aPTT < antiXa”
## # A tibble: 113 × 21
## mrn age_years sex weight height hours_ecmo department_name mode
## <chr> <dbl> <fct> <dbl> <dbl> <dbl> <fct> <fct>
## 1 1000106939 61 Female 63.1 151 111 MIL 5 CCU VA
## 2 1001325965 71 Male 109. 183 90 MIL 5 CTICU VA
## 3 1001491594 51 Male 98.5 178. 189 HRT CARDIAC CARE VA
## 4 1002442790 46 Female 71.5 158. 46 MIL 5 CTICU VV
## 5 1002467250 69 Female 63.5 152. 192 MIL 5 CTICU VA
## 6 1002537539 63 Female 70.4 163. 152 MIL 5 CCU VA
## 7 1002687909 34 Female 84.8 163. 255 MIL 5 CTICU Conver…
## 8 1003328845 72 Female 66 160 179 HRT CARDIAC CARE VA
## 9 1003481283 68 Male 83 168 168 MIL 5 CCU VA
## 10 1003642179 50 Male 128 170 70 HRT CARDIAC CARE VA
## # ℹ 103 more rows
## # ℹ 13 more variables: pump_flow4 <dbl>, pump_flow24 <dbl>,
## # activated_partial_thromboplastin_time <dbl>,
## # heparin_assay_quantitative <dbl>, international_normalization_ratio <dbl>,
## # fibrinogen <dbl>, nearest_bilirubin <dbl>,
## # heparin_infusion_dose_value_units_kg_hr <dbl>, ever_high <lgl>,
## # ever_low <lgl>, ever_disc <lgl>, discord_group <chr>, .groups <chr>| discord_group | n_patients | median_age | iqr_age | male_pct | median_wt | iqr_wt | median_ht | iqr_ht | median_hours | iqr_hours | cticu_pct | ccu_pct | vv_pct | va_pct | median_aptt | median_antixa | median_inr | median_fibrinogen | median_bilirubin | median_hep_dose |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| aPTT < antiXa | 18 | 60 | 19 | 55.56 | 99.05 | 22.5 | 171.6 | 13.30 | 164.5 | 114.75 | 16.67 | 55.56 | 5.56 | 83.33 | 42.90 | 0.36 | 1.2 | 349.0 | 0.7 | 9.5 |
| aPTT > antiXa | 17 | 58 | 16 | 52.94 | 73.00 | 20.4 | 167.6 | 13.00 | 192.0 | 259.00 | 64.71 | 11.76 | 5.88 | 76.47 | 50.55 | 0.17 | 1.2 | 466.5 | 2.0 | 9.0 |
| both types | 1 | 51 | 0 | 100.00 | 98.50 | 0.0 | 177.8 | 0.00 | 189.0 | 0.00 | 0.00 | 0.00 | 0.00 | 100.00 | 61.20 | 0.47 | 1.3 | 210.5 | 0.9 | 11.5 |
| never discordant | 77 | 56 | 21 | 67.53 | 84.80 | 29.7 | 170.2 | 12.93 | 114.0 | 108.00 | 51.95 | 20.78 | 16.88 | 63.64 | 34.00 | 0.15 | 1.1 | 477.5 | 0.8 | 6.0 |
| variable | p_value |
|---|---|
| age_years | 0.8833 |
| weight | 0.0030 |
| height | 0.3919 |
| hours_ecmo | 0.0789 |
| activated_partial_thromboplastin_time | 0.0000 |
| heparin_assay_quantitative | 0.0000 |
| international_normalization_ratio | 0.0974 |
| fibrinogen | 0.1939 |
| nearest_bilirubin | 0.0297 |
| heparin_infusion_dose_value_units_kg_hr | 0.0342 |
| variable | group1 | group2 | p.value |
|---|---|---|---|
| age_years | aPTT > antiXa | aPTT < antiXa | 0.8171 |
| age_years | both types | aPTT < antiXa | 0.7840 |
| age_years | both types | aPTT > antiXa | 0.5627 |
| age_years | never discordant | aPTT < antiXa | 0.7359 |
| age_years | never discordant | aPTT > antiXa | 0.5554 |
| age_years | never discordant | both types | 0.6568 |
| weight | aPTT > antiXa | aPTT < antiXa | 0.0003 |
| weight | both types | aPTT < antiXa | 1.0000 |
| weight | both types | aPTT > antiXa | 0.1770 |
| weight | never discordant | aPTT < antiXa | 0.0169 |
| weight | never discordant | aPTT > antiXa | 0.0248 |
| weight | never discordant | both types | 0.4773 |
| height | aPTT > antiXa | aPTT < antiXa | 0.2033 |
| height | both types | aPTT < antiXa | 0.5822 |
| height | both types | aPTT > antiXa | 0.2888 |
| height | never discordant | aPTT < antiXa | 0.6203 |
| height | never discordant | aPTT > antiXa | 0.2046 |
| height | never discordant | both types | 0.3558 |
| hours_ecmo | aPTT > antiXa | aPTT < antiXa | 0.5747 |
| hours_ecmo | both types | aPTT < antiXa | 0.9474 |
| hours_ecmo | both types | aPTT > antiXa | 1.0000 |
| hours_ecmo | never discordant | aPTT < antiXa | 0.1699 |
| hours_ecmo | never discordant | aPTT > antiXa | 0.0191 |
| hours_ecmo | never discordant | both types | 0.3987 |
| activated_partial_thromboplastin_time | aPTT > antiXa | aPTT < antiXa | 0.1949 |
| activated_partial_thromboplastin_time | both types | aPTT < antiXa | 0.4211 |
| activated_partial_thromboplastin_time | both types | aPTT > antiXa | 0.3333 |
| activated_partial_thromboplastin_time | never discordant | aPTT < antiXa | 0.0008 |
| activated_partial_thromboplastin_time | never discordant | aPTT > antiXa | 0.0000 |
| activated_partial_thromboplastin_time | never discordant | both types | 0.1427 |
| heparin_assay_quantitative | aPTT > antiXa | aPTT < antiXa | 0.0000 |
| heparin_assay_quantitative | both types | aPTT < antiXa | 0.4088 |
| heparin_assay_quantitative | both types | aPTT > antiXa | 0.1209 |
| heparin_assay_quantitative | never discordant | aPTT < antiXa | 0.0000 |
| heparin_assay_quantitative | never discordant | aPTT > antiXa | 0.3838 |
| heparin_assay_quantitative | never discordant | both types | 0.0998 |
| international_normalization_ratio | aPTT > antiXa | aPTT < antiXa | 0.5597 |
| international_normalization_ratio | both types | aPTT < antiXa | 0.2521 |
| international_normalization_ratio | both types | aPTT > antiXa | 0.4346 |
| international_normalization_ratio | never discordant | aPTT < antiXa | 0.1807 |
| international_normalization_ratio | never discordant | aPTT > antiXa | 0.0472 |
| international_normalization_ratio | never discordant | both types | 0.2398 |
| fibrinogen | aPTT > antiXa | aPTT < antiXa | 0.2544 |
| fibrinogen | both types | aPTT < antiXa | 0.9091 |
| fibrinogen | both types | aPTT > antiXa | 0.3333 |
| fibrinogen | never discordant | aPTT < antiXa | 0.1218 |
| fibrinogen | never discordant | aPTT > antiXa | 0.8572 |
| fibrinogen | never discordant | both types | 0.1120 |
| nearest_bilirubin | aPTT > antiXa | aPTT < antiXa | 0.0643 |
| nearest_bilirubin | both types | aPTT < antiXa | 0.9271 |
| nearest_bilirubin | both types | aPTT > antiXa | 0.3333 |
| nearest_bilirubin | never discordant | aPTT < antiXa | 0.6826 |
| nearest_bilirubin | never discordant | aPTT > antiXa | 0.0027 |
| nearest_bilirubin | never discordant | both types | 0.9645 |
| heparin_infusion_dose_value_units_kg_hr | aPTT > antiXa | aPTT < antiXa | 0.3721 |
| heparin_infusion_dose_value_units_kg_hr | both types | aPTT < antiXa | 0.7840 |
| heparin_infusion_dose_value_units_kg_hr | both types | aPTT > antiXa | 0.5623 |
| heparin_infusion_dose_value_units_kg_hr | never discordant | aPTT < antiXa | 0.0071 |
| heparin_infusion_dose_value_units_kg_hr | never discordant | aPTT > antiXa | 0.1874 |
| heparin_infusion_dose_value_units_kg_hr | never discordant | both types | 0.3167 |
Patient-level Demographics, based on FIRST DISCORDANCE TYPE ONLY - descriptive only for consideration in survival analyses (which are based on first discordance type), not planning to report
Note that one patient with both discordance groups is included in this analysis. This patient’s first discordance type was “aPTT < antiXa”
| name | concordant | aPTT > antiXa | aPTT < antiXa |
|---|---|---|---|
| n_patients | 96.00 | 7.00 | 10.00 |
| n_labs | 473.00 | 45.00 | 53.00 |
| median_age | 56.00 | 58.00 | 51.00 |
| iqr_age | 23.00 | 7.00 | 28.00 |
| male_pct | 63.42 | 91.11 | 56.60 |
| median_wt | 86.00 | 83.90 | 98.50 |
| iqr_wt | 30.70 | 13.00 | 9.10 |
| median_ht | 170.20 | 167.60 | 175.00 |
| iqr_ht | 15.20 | 7.70 | 17.80 |
| median_hours | 203.00 | 106.00 | 189.00 |
| iqr_hours | 241.00 | 248.00 | 100.00 |
| cticu_pct | 39.96 | 28.89 | 24.53 |
| ccu_pct | 38.05 | 2.22 | 39.62 |
| hrtccu_pct | 18.82 | 66.67 | 35.85 |
| micu_pct | 2.11 | 0.00 | 0.00 |
| vv_pct | 10.15 | 0.00 | 1.89 |
| va_pct | 67.65 | 100.00 | 96.23 |
| vva_pct | 2.54 | 0.00 | 0.00 |
| convert_pct | 19.66 | 0.00 | 1.89 |
| median_flow4 | 3.10 | 3.20 | 2.80 |
| iqr_flow4 | 0.80 | 0.20 | 1.00 |
| median_flow24 | 3.30 | 3.00 | 3.40 |
| iqr_flow24 | 0.70 | 0.20 | 0.60 |
| median_aptt | 37.40 | 58.10 | 41.80 |
| iqr_aptt | 14.20 | 43.10 | 21.30 |
| median_antixa | 0.18 | 0.30 | 0.38 |
| iqr_antixa | 0.19 | 0.18 | 0.26 |
| median_inr | 1.10 | 1.20 | 1.20 |
| iqr_inr | 0.20 | 0.10 | 0.20 |
| median_fibrinogen | 482.00 | 592.50 | 300.50 |
| iqr_fibrinogen | 349.00 | 246.50 | 153.75 |
| median_heparin_dose | 9.00 | 12.00 | 8.00 |
| iqr_heparin_dose | 6.50 | 6.00 | 4.00 |
| median_nearest_bilirubin | 1.10 | 1.70 | 0.70 |
| iqr_nearest_bilirubin | 1.20 | 1.60 | 0.40 |
Discordance vs Geographic ICU and Support Type - for Connie abstract
##
## concordant aPTT > antiXa aPTT < antiXa
## ccu 271 16 53
## cticu 191 19 5
## micu 10 0 0##
## Fisher's Exact Test for Count Data
##
## data: table_locations
## p-value = 1.245e-06
## alternative hypothesis: two.sided##
## concordant aPTT > antiXa aPTT < antiXa
## Cardiac 365 31 57
## ECPR 32 3 0
## Pulmonary 80 2 1## mode
## unusual_type VA VV
## concordant 338 46
## aPTT > antiXa 30 2
## aPTT < antiXa 48 1##
## Fisher's Exact Test for Count Data
##
## data: .
## p-value = 0.05746
## alternative hypothesis: two.sidedSafety checks to ensure only including labs and complication that were measured between ECMO start and 7 days following ECMO STOP
# safety check to ensure only including labs that were measured after ecmo start or before 7 days following ECMO stop
model_fit_data %>%
filter(lab_result_time < time_on | lab_result_time > complication_7day_assessment)## # A tibble: 0 × 139
## # ℹ 139 variables: study_id <int>, mrn <chr>, patient_name <chr>,
## # time_on <dttm>, time_off <dttm>, time_off_true <dttm>,
## # complication_7day_assessment <dttm>, duration_to_lab <dbl>,
## # patient_id <chr>, run_id <chr>, lab_result_time <dttm>,
## # department_name <fct>, heparin_infusion_dose_value_units_kg_hr <dbl>,
## # concordance_classic <fct>, concordance_low <fct>,
## # concordance_classic_simple <fct>, concordance_low_simple <fct>, …## check shows no data -> all is okay#safety check to ensure only including complications that were measured after ecmo start or before 7 days following ECMO stop
model_fit_data %>%
pivot_longer(
death_date:ich_based_on_imaging_findings,
names_to = "complication",
values_to = "time_complication") %>%
filter(time_complication < time_on | time_complication > complication_7day_assessment)## # A tibble: 0 × 98
## # ℹ 98 variables: study_id <int>, mrn <chr>, patient_name <chr>,
## # time_on <dttm>, time_off <dttm>, time_off_true <dttm>,
## # complication_7day_assessment <dttm>, duration_to_lab <dbl>,
## # patient_id <chr>, run_id <chr>, lab_result_time <dttm>,
## # department_name <fct>, heparin_infusion_dose_value_units_kg_hr <dbl>,
## # concordance_classic <fct>, concordance_low <fct>,
## # concordance_classic_simple <fct>, concordance_low_simple <fct>, …## check shows no data -> all is okayCOMPLICATION ANALYSES
Complications are assessed only if they occur after a lab (discordance) result time and before 7 days after ECMO stop.
Created composite hemorrhagic and thrombotic outcomes, grouped as below.
Complication Composite Groupings
Any Hemorrhagic Complication =
hemorrhage_any
neurologic_intra_extra_parenchymal_cns_hemorrhage_us_or_ct_or_mrihemorrhagic_gi_hemorrhagemajor_bleeding_perfusionhemorrhagic_peripheral_cannulation_site_bleedingpulmonary_pulmonary_hemorrhagehemorrhagic_surgical_site_bleedinghemorrhagic_mediastinal_cannulation_site_bleedinggi_bleeding_present_on_cta_and_confirmed_via_endoscopy***ich_based_on_imaging_findings***
Any Thrombotic Complication =
thrombosis_any
neurologic_cns_infarction_us_or_ct_or_mrilimb_ischemiamechanical_thrombosis_clots_circuit_componentcircuit_thrombosis_perfusiondvt_pe_1_time_and_date***ischemic_stroke_1_time_and_date***acute_limb_ischemia_1_time_and_date***
not included in thrombosis composite as unlikely to be biologically related to discordance:
neurologic_cns_diffuse_ischemia_ct_mrilimb_fasciotomy
*** = complications as reviewed manually by Phil/Connie
Skipped as could violate causal time pathway between discordance -> complication
dvt_pe_2_time_and_date
dvt_pe_3_time_and_date
dvt_pe_4_time_and_date
acute_limb_ischemia_2_time_and_date
acute_limb_ischemia_3_time_and_date
ischemic_stroke_2_time_and_date
ischemic_stroke_3_time_and_date
| complication | aPTT < antiXa | aPTT > antiXa | both types | never discordant |
|---|---|---|---|---|
| any_acute_limb_ischemia_1_time_and_date | 0 | 2 | 0 | 7 |
| any_circuit_thrombosis_perfusion | 3 | 1 | 0 | 15 |
| any_dvt_pe_1_time_and_date | 1 | 2 | 0 | 14 |
| any_gi_bleeding_present_on_cta_and_confirmed_via_endoscopy | 0 | 1 | 0 | 0 |
| any_hemorrhagic_gi_hemorrhage | 0 | 0 | 0 | 1 |
| any_hemorrhagic_mediastinal_cannulation_site_bleeding | 0 | 0 | 0 | 0 |
| any_hemorrhagic_peripheral_cannulation_site_bleeding | 0 | 1 | 0 | 0 |
| any_hemorrhagic_surgical_site_bleeding | 0 | 0 | 0 | 1 |
| any_ich_based_on_imaging_findings | 1 | 2 | 0 | 0 |
| any_ischemic_stroke_1_time_and_date | 0 | 2 | 0 | 2 |
| any_limb_ischemia | 0 | 0 | 0 | 0 |
| any_major_bleeding_perfusion | 0 | 1 | 0 | 1 |
| any_mechanical_thrombosis_clots_circuit_component | 0 | 0 | 0 | 0 |
| any_neurologic_cns_infarction_us_or_ct_or_mri | 0 | 1 | 0 | 0 |
| any_neurologic_intra_extra_parenchymal_cns_hemorrhage_us_or_ct_or_mri | 0 | 0 | 0 | 0 |
| any_pulmonary_pulmonary_hemorrhage | 0 | 0 | 0 | 0 |
| hemorrhage_any | 1 | 4 | 0 | 3 |
| thrombosis_any | 4 | 6 | 0 | 34 |
patient_model_data <- patient_level_complications %>%
left_join(patient_demographics, by = "mrn") %>%
filter(discord_group != "both types") %>%
mutate(
discord_group = fct_relevel(discord_group, "never discordant"),
icu_type = case_when(
department_name == "MIL 5 CCU" ~ "ccu",
department_name == "MIL 5 CTICU" ~ "cticu",
department_name == "MIL 4 SICU" ~ "other",
department_name == "MIL OPERATING ROOM" ~ "other",
department_name == "HRT CARDIAC CARE" ~ "ccu",
department_name == "MIL 4 MICU A" ~ "micu",
department_name == "MIL CARDIAC CATH" ~ "other",
department_name == "PBY ADULT EMERGENCY" ~ "other",
department_name == "MIL 4 MICU B" ~ "micu",
department_name == "GBG 4 W CT ICU" ~ "other",
department_name == "MIL 9 HUDSON" ~ "other",
department_name == "POST-DISCHARGE" ~ "other",
TRUE ~ NA)) %>%
filter(icu_type != "other")0. patient-level analyses
fibrinogen only measured in 38 patients in this model
nearest_bilirubin, international_normalization_ratio, weight, hours_ecmohemorrhage_patient_glm <- glm(
hemorrhage_any ~ discord_group,
data = patient_model_data,
family = binomial)
hemorrhage_patient_glm %>% broom::tidy()## # A tibble: 3 × 5
## term estimate std.error statistic p.value
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) -3.16 0.589 -5.37 0.0000000796
## 2 discord_groupaPTT < antiXa 0.331 1.19 0.279 0.780
## 3 discord_groupaPTT > antiXa 2.07 0.825 2.50 0.0123thrombosis_patient_glm <- glm(
thrombosis_any ~ discord_group,
data = patient_model_data,
family = binomial)
thrombosis_patient_glm %>% broom::tidy(exponentiate = TRUE, conf.int = TRUE)## # A tibble: 3 × 7
## term estimate std.error statistic p.value conf.low conf.high
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 0.805 0.234 -0.928 0.353 0.506 1.27
## 2 discord_groupaPTT < a… 0.355 0.613 -1.69 0.0913 0.0937 1.10
## 3 discord_groupaPTT > a… 0.745 0.567 -0.518 0.604 0.233 2.22hemorrhage_group_sick = glm(hemorrhage_any ~ discord_group + pump_flow24 + heparin_infusion_dose_value_units_kg_hr, data = patient_model_data, family = binomial)
hemorrhage_group_sick %>% broom::tidy(exponentiate = TRUE, conf.int = TRUE)## # A tibble: 5 × 7
## term estimate std.error statistic p.value conf.low conf.high
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 0.000245 3.46 -2.40 0.0163 9.55e-8 0.110
## 2 discord_groupaPTT < a… 0.935 1.28 -0.0530 0.958 3.93e-2 9.14
## 3 discord_groupaPTT > a… 8.04 0.867 2.40 0.0162 1.47e+0 49.3
## 4 pump_flow24 3.73 0.966 1.36 0.173 6.09e-1 29.8
## 5 heparin_infusion_dose… 1.11 0.101 1.03 0.303 9.13e-1 1.37thrombosis_group_sick = glm(thrombosis_any ~ discord_group + pump_flow24 + heparin_infusion_dose_value_units_kg_hr, data = patient_model_data, family = binomial)
thrombosis_group_sick %>% broom::tidy(exponentiate = TRUE, conf.int = TRUE)## # A tibble: 5 × 7
## term estimate std.error statistic p.value conf.low conf.high
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 0.363 1.44 -0.705 0.481 0.0206 6.00
## 2 discord_groupaPTT < a… 0.331 0.713 -1.55 0.121 0.0683 1.22
## 3 discord_groupaPTT > a… 0.727 0.615 -0.519 0.604 0.203 2.37
## 4 pump_flow24 1.41 0.444 0.771 0.441 0.592 3.43
## 5 heparin_infusion_dose… 0.960 0.0550 -0.736 0.462 0.860 1.07review
discord_grp_glm = glm(discord_group ~ international_normalization_ratio + weight + heparin_infusion_dose_value_units_kg_hr + hours_ecmo, data = patient_model_data, family = binomial)
discord_grp_glm##
## Call: glm(formula = discord_group ~ international_normalization_ratio +
## weight + heparin_infusion_dose_value_units_kg_hr + hours_ecmo,
## family = binomial, data = patient_model_data)
##
## Coefficients:
## (Intercept)
## -3.540743
## international_normalization_ratio
## 0.960743
## weight
## -0.001890
## heparin_infusion_dose_value_units_kg_hr
## 0.161465
## hours_ecmo
## 0.001955
##
## Degrees of Freedom: 98 Total (i.e. Null); 94 Residual
## (9 observations deleted due to missingness)
## Null Deviance: 123.1
## Residual Deviance: 109.5 AIC: 119.5discord_grp_glm %>% broom::tidy()## # A tibble: 5 × 5
## term estimate std.error statistic p.value
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) -3.54 1.56 -2.26 0.0237
## 2 international_normalization_ratio 0.961 1.01 0.954 0.340
## 3 weight -0.00189 0.00924 -0.205 0.838
## 4 heparin_infusion_dose_value_units_kg_hr 0.161 0.0580 2.78 0.00541
## 5 hours_ecmo 0.00195 0.00134 1.46 0.1431. unadjusted analysis - will not report
## # A tibble: 3 × 4
## unusual_type n hemorrhage_rate thrombosis_rate
## <fct> <int> <dbl> <dbl>
## 1 concordant 477 7.13 29.8
## 2 aPTT > antiXa 36 11.1 25
## 3 aPTT < antiXa 58 1.72 24.1##
## Call: glm(formula = hemorrhage_any ~ unusual_type, family = binomial,
## data = model_fit_complications)
##
## Coefficients:
## (Intercept) unusual_typeaPTT > antiXa
## -2.5672 0.4878
## unusual_typeaPTT < antiXa
## -1.4758
##
## Degrees of Freedom: 570 Total (i.e. Null); 568 Residual
## Null Deviance: 284.6
## Residual Deviance: 280.3 AIC: 286.3##
## Call:
## glm(formula = hemorrhage_any ~ unusual_type, family = binomial,
## data = model_fit_complications)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.4854 -0.3846 -0.3846 -0.3846 2.8497
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -2.5672 0.1780 -14.426 <2e-16 ***
## unusual_typeaPTT > antiXa 0.4878 0.5594 0.872 0.383
## unusual_typeaPTT < antiXa -1.4758 1.0241 -1.441 0.150
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 284.61 on 570 degrees of freedom
## Residual deviance: 280.33 on 568 degrees of freedom
## AIC: 286.33
##
## Number of Fisher Scoring iterations: 6##
## Call: glm(formula = thrombosis_any ~ unusual_type, family = binomial,
## data = model_fit_complications)
##
## Coefficients:
## (Intercept) unusual_typeaPTT > antiXa
## -0.8583 -0.2403
## unusual_typeaPTT < antiXa
## -0.2868
##
## Degrees of Freedom: 570 Total (i.e. Null); 568 Residual
## Null Deviance: 686.6
## Residual Deviance: 685.5 AIC: 691.5##
## Call:
## glm(formula = thrombosis_any ~ unusual_type, family = binomial,
## data = model_fit_complications)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.8407 -0.8407 -0.8407 1.5567 1.6861
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.8583 0.1001 -8.571 <2e-16 ***
## unusual_typeaPTT > antiXa -0.2403 0.3977 -0.604 0.546
## unusual_typeaPTT < antiXa -0.2868 0.3228 -0.889 0.374
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 686.60 on 570 degrees of freedom
## Residual deviance: 685.49 on 568 degrees of freedom
## AIC: 691.49
##
## Number of Fisher Scoring iterations: 42. adjusted analysis = PRIMARY OUTCOME RESULT
Variables which would be important to adjust for based on opinion / intentional selection:
duration_to_lab(time_on -> lab_result_time)international_normalization_ratiofibrinogennearest_bilirubinweight- significant as identified with between-group testingheparin_infusion_dose_value_units_kg_hr
Variables which specifically were not incorporated into the model as they were not between-group significant or had other good reasons:
pump_flow4andpump_flow24- specifically not chosen as they were not significantly different between groups in descriptive analyses, despite potentially functioning as marker of baseline disease severity (higher flow, sicker). We do not know illness severity and underlying disease process in this retrospective cohort study.mode-> inclusion of this model caused model to possibly overfit AND our underlying premise is that discordance matters, regardless of why you are on ECMO or how it is cannulated
##
## Call:
## glm(formula = hemorrhage_any ~ unusual_type + duration_to_lab +
## international_normalization_ratio + fibrinogen + nearest_bilirubin +
## weight + heparin_infusion_dose_value_units_kg_hr, family = binomial,
## data = model_fit_complications)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.7662 -0.3380 -0.2056 -0.1233 3.2015
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 5.437e-01 3.189e+00 0.170 0.86462
## unusual_typeaPTT > antiXa -1.629e+01 1.751e+03 -0.009 0.99258
## unusual_typeaPTT < antiXa -4.178e-01 1.301e+00 -0.321 0.74816
## duration_to_lab -1.206e-03 4.098e-03 -0.294 0.76857
## international_normalization_ratio 1.358e+00 2.296e+00 0.592 0.55404
## fibrinogen 2.194e-05 1.954e-03 0.011 0.99104
## nearest_bilirubin -4.498e-01 4.408e-01 -1.020 0.30753
## weight -6.816e-02 2.174e-02 -3.135 0.00172
## heparin_infusion_dose_value_units_kg_hr 1.732e-01 9.296e-02 1.863 0.06248
##
## (Intercept)
## unusual_typeaPTT > antiXa
## unusual_typeaPTT < antiXa
## duration_to_lab
## international_normalization_ratio
## fibrinogen
## nearest_bilirubin
## weight **
## heparin_infusion_dose_value_units_kg_hr .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 86.459 on 167 degrees of freedom
## Residual deviance: 61.610 on 159 degrees of freedom
## (403 observations deleted due to missingness)
## AIC: 79.61
##
## Number of Fisher Scoring iterations: 16| GVIF | Df | GVIF^(1/(2*Df)) | |
|---|---|---|---|
| unusual_type | 1.355896 | 2 | 1.079087 |
| duration_to_lab | 1.755864 | 1 | 1.325090 |
| international_normalization_ratio | 1.304064 | 1 | 1.141956 |
| fibrinogen | 1.407007 | 1 | 1.186173 |
| nearest_bilirubin | 1.147931 | 1 | 1.071416 |
| weight | 1.595967 | 1 | 1.263316 |
| heparin_infusion_dose_value_units_kg_hr | 1.836584 | 1 | 1.355206 |
## # A tibble: 9 × 8
## term OR conf.low conf.high p_raw p_fdr p_holm p_bonf
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 1.72e+0 0.00588 2.22e+ 3 0.865 0.993 1 1
## 2 unusual_typeaPTT > an… 8.44e-8 NA 7.69e+56 0.993 0.993 1 1
## 3 unusual_typeaPTT < an… 6.59e-1 0.0276 7.14e+ 0 0.748 0.993 1 1
## 4 duration_to_lab 9.99e-1 0.990 1.01e+ 0 0.769 0.993 1 1
## 5 international_normali… 3.89e+0 0.0134 1.70e+ 2 0.554 0.993 1 1
## 6 fibrinogen 1.00e+0 0.996 1.00e+ 0 0.991 0.993 1 1
## 7 nearest_bilirubin 6.38e-1 0.209 1.21e+ 0 0.308 0.923 1 1
## 8 weight 9.34e-1 0.890 9.71e- 1 0.00172 0.0155 0.0155 0.0155
## 9 heparin_infusion_dose… 1.19e+0 1.00 1.45e+ 0 0.0625 0.281 0.500 0.562## [1] TRUE##
## FALSE TRUE
## concordant 443 34
## aPTT > antiXa 32 4
## aPTT < antiXa 57 1
| term | estimate | OR | conf.low | conf.high | p_raw | p_fdr | p_holm | p_bonf |
|---|---|---|---|---|---|---|---|---|
| (Intercept) | -0.4517335 | 0.6365238 | 0.0055107 | 394.8943525 | 0.8573556 | 0.9645251 | 1.0000000 | 1.0000000 |
| unusual_typeaPTT > antiXa | -1.1943994 | 0.3028858 | 0.0019908 | 4.7298196 | 0.4328422 | 0.8319527 | 1.0000000 | 1.0000000 |
| unusual_typeaPTT < antiXa | -0.2020069 | 0.8170893 | 0.0680678 | 6.1958257 | 0.8512867 | 0.9645251 | 1.0000000 | 1.0000000 |
| duration_to_lab | -0.0008927 | 0.9991077 | 0.9914045 | 1.0066575 | 0.8082157 | 0.9645251 | 1.0000000 | 1.0000000 |
| international_normalization_ratio | 1.7131475 | 5.5463913 | 0.0361873 | 90.6128738 | 0.3617986 | 0.8319527 | 1.0000000 | 1.0000000 |
| fibrinogen | -0.0000540 | 0.9999460 | 0.9962827 | 1.0036019 | 0.9762417 | 0.9762417 | 1.0000000 | 1.0000000 |
| nearest_bilirubin | -0.2232760 | 0.7998941 | 0.3003558 | 1.2897960 | 0.4621959 | 0.8319527 | 1.0000000 | 1.0000000 |
| weight | -0.0575024 | 0.9441196 | 0.9050462 | 0.9775113 | 0.0009199 | 0.0082792 | 0.0082792 | 0.0082792 |
| heparin_infusion_dose_value_units_kg_hr | 0.1432434 | 1.1540107 | 0.9887063 | 1.3823250 | 0.0701087 | 0.3154892 | 0.5608697 | 0.6309784 |
##
## Call:
## glm(formula = thrombosis_any ~ unusual_type + international_normalization_ratio +
## heparin_infusion_dose_value_units_kg_hr + duration_to_lab +
## fibrinogen + nearest_bilirubin + weight, family = binomial,
## data = model_fit_complications)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -2.2317 -0.8228 -0.4232 0.9428 2.4156
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 4.201e-01 1.858e+00 0.226 0.821126
## unusual_typeaPTT > antiXa -1.488e+01 1.056e+03 -0.014 0.988756
## unusual_typeaPTT < antiXa 1.348e+00 5.398e-01 2.497 0.012509
## international_normalization_ratio -5.654e-01 1.228e+00 -0.460 0.645206
## heparin_infusion_dose_value_units_kg_hr -2.186e-01 6.284e-02 -3.478 0.000505
## duration_to_lab -7.959e-03 3.629e-03 -2.193 0.028276
## fibrinogen -8.438e-04 9.307e-04 -0.907 0.364581
## nearest_bilirubin 2.562e-02 1.347e-01 0.190 0.849085
## weight 2.474e-02 1.010e-02 2.449 0.014327
##
## (Intercept)
## unusual_typeaPTT > antiXa
## unusual_typeaPTT < antiXa *
## international_normalization_ratio
## heparin_infusion_dose_value_units_kg_hr ***
## duration_to_lab *
## fibrinogen
## nearest_bilirubin
## weight *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 202.82 on 167 degrees of freedom
## Residual deviance: 162.36 on 159 degrees of freedom
## (403 observations deleted due to missingness)
## AIC: 180.36
##
## Number of Fisher Scoring iterations: 15## # A tibble: 9 × 7
## term estimate std.error statistic p.value conf.low conf.high
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 1.52e+0 1.86e+0 0.226 8.21e-1 0.0413 7.45e+ 1
## 2 unusual_typeaPTT > an… 3.44e-7 1.06e+3 -0.0141 9.89e-1 NA 3.54e+40
## 3 unusual_typeaPTT < an… 3.85e+0 5.40e-1 2.50 1.25e-2 1.37 1.16e+ 1
## 4 international_normali… 5.68e-1 1.23e+0 -0.460 6.45e-1 0.0384 6.00e+ 0
## 5 heparin_infusion_dose… 8.04e-1 6.28e-2 -3.48 5.05e-4 0.705 9.03e- 1
## 6 duration_to_lab 9.92e-1 3.63e-3 -2.19 2.83e-2 0.985 9.99e- 1
## 7 fibrinogen 9.99e-1 9.31e-4 -0.907 3.65e-1 0.997 1.00e+ 0
## 8 nearest_bilirubin 1.03e+0 1.35e-1 0.190 8.49e-1 0.769 1.32e+ 0
## 9 weight 1.03e+0 1.01e-2 2.45 1.43e-2 1.01 1.05e+ 0| GVIF | Df | GVIF^(1/(2*Df)) | |
|---|---|---|---|
| unusual_type | 1.258849 | 2 | 1.059238 |
| international_normalization_ratio | 1.408611 | 1 | 1.186849 |
| heparin_infusion_dose_value_units_kg_hr | 1.503484 | 1 | 1.226166 |
| duration_to_lab | 1.323139 | 1 | 1.150278 |
| fibrinogen | 1.322245 | 1 | 1.149889 |
| nearest_bilirubin | 1.418127 | 1 | 1.190851 |
| weight | 1.346531 | 1 | 1.160401 |
## # A tibble: 9 × 8
## term OR conf.low conf.high p_raw p_fdr p_holm p_bonf
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 1.52e+0 0.0413 7.45e+ 1 8.21e-1 0.955 1 1
## 2 unusual_typeaPTT >… 3.44e-7 NA 3.54e+40 9.89e-1 0.989 1 1
## 3 unusual_typeaPTT <… 3.85e+0 1.37 1.16e+ 1 1.25e-2 0.0430 0.100 0.113
## 4 international_norm… 5.68e-1 0.0384 6.00e+ 0 6.45e-1 0.955 1 1
## 5 heparin_infusion_d… 8.04e-1 0.705 9.03e- 1 5.05e-4 0.00454 0.00454 0.00454
## 6 duration_to_lab 9.92e-1 0.985 9.99e- 1 2.83e-2 0.0636 0.170 0.254
## 7 fibrinogen 9.99e-1 0.997 1.00e+ 0 3.65e-1 0.656 1 1
## 8 nearest_bilirubin 1.03e+0 0.769 1.32e+ 0 8.49e-1 0.955 1 1
## 9 weight 1.03e+0 1.01 1.05e+ 0 1.43e-2 0.0430 0.100 0.129## [1] TRUE##
## FALSE TRUE
## concordant 335 142
## aPTT > antiXa 27 9
## aPTT < antiXa 44 14
| term | estimate | OR | conf.low | conf.high | p_raw | p_fdr | p_holm | p_bonf |
|---|---|---|---|---|---|---|---|---|
| (Intercept) | 0.3112238 | 1.3650946 | 0.0489856 | 50.5729181 | 0.8542002 | 0.8542002 | 1.0000000 | 1.0000000 |
| unusual_typeaPTT > antiXa | -0.8812509 | 0.4142644 | 0.0029282 | 5.0395611 | 0.5416406 | 0.8124609 | 1.0000000 | 1.0000000 |
| unusual_typeaPTT < antiXa | 1.2449190 | 3.4726536 | 1.2893237 | 9.9483446 | 0.0136282 | 0.0408845 | 0.0953973 | 0.1226536 |
| duration_to_lab | -0.0073757 | 0.9926514 | 0.9854777 | 0.9991030 | 0.0247555 | 0.0556998 | 0.1485328 | 0.2227992 |
| international_normalization_ratio | -0.4233769 | 0.6548318 | 0.0549153 | 5.3981260 | 0.6915597 | 0.8542002 | 1.0000000 | 1.0000000 |
| fibrinogen | -0.0008607 | 0.9991397 | 0.9973761 | 1.0009172 | 0.3391874 | 0.6105373 | 1.0000000 | 1.0000000 |
| nearest_bilirubin | 0.0386472 | 1.0394037 | 0.7886354 | 1.3062943 | 0.7627375 | 0.8542002 | 1.0000000 | 1.0000000 |
| weight | 0.0223096 | 1.0225603 | 1.0054031 | 1.0444797 | 0.0081549 | 0.0366972 | 0.0652394 | 0.0733944 |
| heparin_infusion_dose_value_units_kg_hr | -0.1996967 | 0.8189791 | 0.7221164 | 0.9152349 | 0.0002459 | 0.0022128 | 0.0022128 | 0.0022128 |
Sensitivity Analysis in which Discordance is Grouped - not to report
hemorrhage_simple = glm(hemorrhage_any ~ discordance_any + duration_to_lab + nearest_bilirubin + fibrinogen + international_normalization_ratio + weight + heparin_infusion_dose_value_units_kg_hr, data = model_fit_complications, family = binomial)
hemorrhage_simple %>% broom::tidy(exponentiate = TRUE, conf.int = TRUE)## # A tibble: 8 × 7
## term estimate std.error statistic p.value conf.low conf.high
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 1.44 3.15 0.116 0.908 0.00507 1826.
## 2 discordance_anydiscor… 0.351 1.25 -0.836 0.403 0.0154 3.27
## 3 duration_to_lab 1.00 0.00402 -0.102 0.919 0.991 1.01
## 4 nearest_bilirubin 0.635 0.451 -1.01 0.314 0.204 1.21
## 5 fibrinogen 1.00 0.00191 -0.152 0.879 0.996 1.00
## 6 international_normali… 4.12 2.25 0.629 0.529 0.0138 173.
## 7 weight 0.938 0.0209 -3.07 0.00218 0.896 0.974
## 8 heparin_infusion_dose… 1.17 0.0900 1.74 0.0818 0.989 1.41thrombosis_simple = glm(thrombosis_any ~ discordance_any + duration_to_lab + international_normalization_ratio + weight + heparin_infusion_dose_value_units_kg_hr + nearest_bilirubin + fibrinogen, data = model_fit_complications, family = binomial)
thrombosis_simple %>% broom::tidy(exponentiate = TRUE, conf.int = TRUE)## # A tibble: 8 × 7
## term estimate std.error statistic p.value conf.low conf.high
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 0.994 1.86 -0.00297 9.98e-1 0.0268 48.6
## 2 discordance_anydiscor… 2.96 0.513 2.12 3.42e-2 1.10 8.34
## 3 duration_to_lab 0.993 0.00349 -1.90 5.77e-2 0.986 1.00
## 4 international_normali… 0.591 1.23 -0.428 6.69e-1 0.0397 6.26
## 5 weight 1.03 0.0103 2.87 4.17e-3 1.01 1.05
## 6 heparin_infusion_dose… 0.787 0.0629 -3.82 1.36e-4 0.690 0.884
## 7 nearest_bilirubin 1.05 0.133 0.367 7.13e-1 0.791 1.34
## 8 fibrinogen 0.999 0.000927 -1.04 2.98e-1 0.997 1.00Sensitivity Analysis in which outcomes are grouped
complication_simple = glm(coag_event_any ~ unusual_type + duration_to_lab + nearest_bilirubin + fibrinogen + international_normalization_ratio + weight + heparin_infusion_dose_value_units_kg_hr, data = model_fit_complications, family = binomial)
complication_simple %>% broom::tidy()## # A tibble: 8 × 5
## term estimate std.error statistic p.value
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) -2.66e+ 1 545666. -4.87e- 5 1.00
## 2 unusual_typeaPTT < antiXa 8.01e-15 124884. 6.41e-20 1
## 3 duration_to_lab -4.68e-17 625. -7.50e-20 1
## 4 nearest_bilirubin -2.90e-15 50094. -5.79e-20 1
## 5 fibrinogen -8.97e-18 298. -3.01e-20 1
## 6 international_normalization_ratio 2.83e-14 349714. 8.09e-20 1
## 7 weight -1.13e-16 2719. -4.15e-20 1
## 8 heparin_infusion_dose_value_units_kg_hr 1.42e-16 11855. 1.20e-20 1Standard logistic regression had high separation based on discordance type, likely due to uncommon outcome among a relatively small sample size with several critical covariables. Thus, Firth penalized logistic regression used to compute odds ratios and p-values.
3. Survival Analyses
- Note that all patients retained in this model. The patient (MRN 1001491594) with multiple types of discordance is therefore included but will be analyzed by whatever happened first. Specifically, this patient’s first instance of discordance was as follows: “2024-08-20 04:09:00 aPTT < antiXa”, thus is included in the aPTT < antiXa survival analysis.
| first_discord_type | n |
|---|---|
| concordant | 77 |
| aPTT > antiXa | 17 |
| aPTT < antiXa | 19 |
| event_throm | n |
|---|---|
| 0 | 69 |
| 1 | 44 |
| event_bleed | n |
|---|---|
| 0 | 106 |
| 1 | 7 |
| discord_group | n |
|---|---|
| concordant | 113 |
| discord_group | n |
|---|---|
| concordant | 113 |
| mrn | time_on | first_discord_time | first_discord_type | clot_first | censor_time |
|---|---|---|---|---|---|
| 1201199342 | 2024-09-09 17:56:00 | 2024-09-26 12:28:00 | aPTT > antiXa | 2024-09-23 22:16:00 | 2024-10-12 17:50:00 |
| mrn | time_on | first_discord_time | first_discord_type | bleed_first | censor_time |
|---|
| mrn | time_on | censor_time | first_discord_time | first_discord_type |
|---|
| mrn | discord_group | switch_time |
|---|---|---|
| 1101397979 | aPTT < antiXa | 0.0125000 |
| 1100547209 | aPTT > antiXa | 0.0854167 |
| 1401591974 | aPTT < antiXa | 0.1076389 |
| 1011043983 | aPTT > antiXa | 0.3354167 |
| 1102216636 | aPTT < antiXa | 0.4888889 |
| 1009544965 | aPTT < antiXa | 0.5597222 |
| 1401327774 | aPTT > antiXa | 0.8465278 |
| 1010147999 | aPTT > antiXa | 1.1416667 |
| 1401183967 | aPTT > antiXa | 1.2652778 |
| 1001491594 | aPTT < antiXa | 1.3305556 |
| 1003642179 | aPTT < antiXa | 1.4020833 |
| 1400032239 | aPTT > antiXa | 1.5722222 |
| 1401696145 | aPTT < antiXa | 1.5930556 |
| 1401639110 | aPTT > antiXa | 1.6680556 |
| 1011077110 | aPTT > antiXa | 1.8166667 |
| 1201555171 | aPTT < antiXa | 2.2041667 |
| 1401689750 | aPTT < antiXa | 2.3638889 |
| 1006905710 | aPTT > antiXa | 2.6826389 |
| 1401697961 | aPTT < antiXa | 3.0881944 |
| 1002467250 | aPTT > antiXa | 3.3597222 |
| 1401682296 | aPTT < antiXa | 3.8819444 |
| 1009477257 | aPTT > antiXa | 4.2250000 |
| 1010696012 | aPTT < antiXa | 4.2479167 |
| 1401611898 | aPTT < antiXa | 4.3090278 |
| 1400969240 | aPTT > antiXa | 4.3513889 |
| 1401468984 | aPTT > antiXa | 5.1944444 |
| 1401554006 | aPTT < antiXa | 5.4881944 |
| 1401707171 | aPTT < antiXa | 6.4451389 |
| 1401991907 | aPTT < antiXa | 6.4826389 |
| 1401848048 | aPTT < antiXa | 7.7048611 |
| 1401766288 | aPTT > antiXa | 7.8972222 |
| 1401796928 | aPTT < antiXa | 9.1236111 |
| 1007129753 | aPTT > antiXa | 13.2458333 |
| 1401683148 | aPTT < antiXa | 20.9833333 |
| 1102078503 | aPTT > antiXa | 26.5215278 |
## Call: survfit(formula = Surv(start, stop, event) ~ discord_group, data = tv_thrombosis)
##
## records n.max n.start events median 0.95LCL 0.95UCL
## discord_group=concordant 113 113 113 35 14.6 13.62 NA
## discord_group=aPTT > antiXa 16 12 8 5 22.2 8.21 NA
## discord_group=aPTT < antiXa 19 15 10 4 NA NA NANote that the static curve risk table counts align perfectly with counts of patient by development of any discordance in the above Demographic descriptive comparisons.
Note also that the 113 patients all start with concordant values. Some patients develop discordance later, and then they start 0 time recategorizes them as “0”. The time-varying curves are aligned by interval start times, not global time 0. So, among patients with aPTT < antiXa in the thrombus curve, n.start indicates that 10 patients are “at risk” at the earliest time for this group, based on days since ECMO cannulation. (Though their start time is based on development of discordance, plotted from interval start time, and in risk table based on days since ECMO cannulation).
Note also that this approach is only testing the impact of first discordance. Repeated episodes of discordance, or discordance severity, not being analyzed. Thus, direct comparison with logistic regression is challenging.
Time varying Cox model - SECONDARY OUTCOME RESULT
| term | std.error | robust.se | statistic | p.value | conf.low | conf.high | HR |
|---|---|---|---|---|---|---|---|
| discord_groupaPTT > antiXa | 1.1026 | 1.0399 | 0.3890 | 0.6973 | 0.1952 | 11.5043 | 1.4985 |
| discord_groupaPTT < antiXa | 0.7093 | 0.7447 | 1.1105 | 0.2668 | 0.5312 | 9.8399 | 2.2863 |
| international_normalization_ratio | 1.4593 | 1.1828 | -0.1231 | 0.9020 | 0.0851 | 8.7820 | 0.8645 |
| fibrinogen | 0.0016 | 0.0016 | 0.3960 | 0.6921 | 0.9976 | 1.0037 | 1.0006 |
| heparin_infusion_dose_value_units_kg_hr | 0.0842 | 0.0923 | -1.4960 | 0.1347 | 0.7268 | 1.0438 | 0.8710 |
| nearest_bilirubin | 0.2219 | 0.2174 | 0.0963 | 0.9233 | 0.6669 | 1.5637 | 1.0212 |
| weight | 0.0096 | 0.0079 | -0.8052 | 0.4207 | 0.9784 | 1.0091 | 0.9937 |
| term | std.error | robust.se | statistic | p.value | conf.low | conf.high | HR |
|---|---|---|---|---|---|---|---|
| discord_groupaPTT > antiXa | 2.1887 | 1.1899 | -0.3587 | 0.7198 | 0.0633 | 6.7221 | 0.6526 |
| discord_groupaPTT < antiXa | 1.6402 | 1.4396 | 0.9045 | 0.3657 | 0.2188 | 61.7883 | 3.6772 |
| international_normalization_ratio | 2.9467 | 2.5123 | 0.7168 | 0.4735 | 0.0440 | 832.9776 | 6.0553 |
| fibrinogen | 0.0037 | 0.0029 | 0.8209 | 0.4117 | 0.9967 | 1.0081 | 1.0024 |
| nearest_bilirubin | 0.8974 | 0.6407 | -1.1718 | 0.2413 | 0.1345 | 1.6569 | 0.4720 |
| weight | 0.0320 | 0.0170 | -2.2556 | 0.0241 | 0.9310 | 0.9950 | 0.9625 |
| heparin_infusion_dose_value_units_kg_hr | 0.1882 | 0.2282 | 0.3671 | 0.7136 | 0.6953 | 1.7006 | 1.0874 |
4. Tertiary outcomes, i.e. those which may not be solely thrombotic or hemorrhagic, but might be associated with discordance - for publication
death_datemechanical_circuit_changerenal_creatinine_3_0renal_creatinine_1_5_3_0renal_renal_replacement_therapy_requiredmetabolic_moderate_hemolysismajor_hemolysis_perfusion
Note: composited renal failure types into
renal_failure_any
Ultimately there are very few of these outcomes (except renal failure and death).
## # A tibble: 6 × 2
## # Groups: complication [6]
## complication n
## <chr> <int>
## 1 death_date 4
## 2 major_hemolysis_perfusion 1
## 3 mechanical_circuit_change 3
## 4 renal_creatinine_1_5_3_0 4
## 5 renal_creatinine_3_0 6
## 6 renal_renal_replacement_therapy_required 3## # A tibble: 3 × 8
## term OR conf.low conf.high p_raw p_fdr p_holm p_bonf
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 1.12e-1 8.20e- 2 1.49e-1 5.91e-47 1.77e-46 1.77e-46 1.77e-46
## 2 unusual_typea… 1.12e+0 3.23e- 1 2.97e+0 8.41e- 1 9.85e- 1 1 e+ 0 1 e+ 0
## 3 unusual_typea… 7.73e-8 1.71e-131 3.81e+8 9.85e- 1 9.85e- 1 1 e+ 0 1 e+ 0| term | OR | conf.low | conf.high | p_raw | p_fdr | p_holm | p_bonf |
|---|---|---|---|---|---|---|---|
| (Intercept) | 2.990310e+22 | 2.127065e+09 | 3.222093e+53 | 0.0215303 | 0.0968863 | 0.1924165 | 0.1937726 |
| unusual_typeaPTT > antiXa | 7.999556e+07 | NA | Inf | 0.9993079 | 0.9993079 | 1.0000000 | 1.0000000 |
| unusual_typeaPTT < antiXa | 1.404830e-02 | NA | 1.359008e+195 | 0.9992693 | 0.9993079 | 1.0000000 | 1.0000000 |
| duration_to_lab | 1.039519e+00 | 1.005664e+00 | 1.131491e+00 | 0.1438864 | 0.2158296 | 0.5755456 | 1.0000000 |
| international_normalization_ratio | 0.000000e+00 | 0.000000e+00 | 1.790000e-05 | 0.0213796 | 0.0968863 | 0.1924165 | 0.1924165 |
| fibrinogen | 1.007144e+00 | 9.901896e-01 | 1.033750e+00 | 0.4558495 | 0.5860922 | 1.0000000 | 1.0000000 |
| nearest_bilirubin | 2.191107e-01 | 1.393990e-02 | 9.899436e-01 | 0.1129538 | 0.2033169 | 0.5647691 | 1.0000000 |
| weight | 8.528896e-01 | 6.663353e-01 | 9.560630e-01 | 0.0480300 | 0.1440901 | 0.3362103 | 0.4322704 |
| heparin_infusion_dose_value_units_kg_hr | 3.029280e-02 | 3.350000e-05 | 2.757105e-01 | 0.0817480 | 0.1839331 | 0.4904882 | 0.7357323 |
##
## Call:
## glm(formula = renal_failure_any ~ unusual_type, family = binomial,
## data = renal_failure_df)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.3144 -0.3144 -0.3144 -0.3144 2.6771
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -2.9826 0.2137 -13.955 <2e-16 ***
## unusual_typeaPTT > antiXa -0.5727 1.0365 -0.553 0.581
## unusual_typeaPTT < antiXa -15.5835 856.4644 -0.018 0.985
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 199.10 on 570 degrees of freedom
## Residual deviance: 193.48 on 568 degrees of freedom
## AIC: 199.48
##
## Number of Fisher Scoring iterations: 17| term | OR | conf.low | conf.high | p_raw | p_fdr | p_holm | p_bonf |
|---|---|---|---|---|---|---|---|
| (Intercept) | 0.0881927 | 0.0001785 | 2.519508e+02 | 0.4498842 | 0.6748263 | 1.0000000 | 1.0000000 |
| unusual_typeaPTT > antiXa | 0.0000000 | NA | 7.101530e+298 | 0.9981573 | 0.9981573 | 1.0000000 | 1.0000000 |
| unusual_typeaPTT < antiXa | 0.0000001 | NA | 3.038054e+108 | 0.9953682 | 0.9981573 | 1.0000000 | 1.0000000 |
| duration_to_lab | 0.9613701 | 0.9172843 | 9.908044e-01 | 0.0414353 | 0.3729176 | 0.3729176 | 0.3729176 |
| international_normalization_ratio | 1.9253381 | 0.0046194 | 1.898370e+02 | 0.7767452 | 0.9981573 | 1.0000000 | 1.0000000 |
| fibrinogen | 1.0023813 | 0.9977075 | 1.007400e+00 | 0.3205538 | 0.6748263 | 1.0000000 | 1.0000000 |
| nearest_bilirubin | 0.6416235 | 0.1253764 | 1.446001e+00 | 0.4309281 | 0.6748263 | 1.0000000 | 1.0000000 |
| weight | 1.0296237 | 0.9824900 | 1.078238e+00 | 0.1776735 | 0.5330206 | 1.0000000 | 1.0000000 |
| heparin_infusion_dose_value_units_kg_hr | 0.7607272 | 0.5090479 | 1.067893e+00 | 0.1365768 | 0.5330206 | 1.0000000 | 1.0000000 |
5. Do counts of discordance differ when using classical (range-based) discordance categorization versus those identified with polynomial model? - CONSIDER FOR ADDITIONAL ABSTRACT
Note that this analysis eliminated 1 lab results in which aPTT was above typical goal rangs and time-matched anti-Xa was below typical goal range and 1 lab result in which aPTT was below typical goal range and time-matched anti-Xa was above typical goal range.
## # A tibble: 5 × 2
## # Groups: concordance_classic_simple [5]
## concordance_classic_simple n
## <fct> <int>
## 1 concordant 379
## 2 discordant aPTT below antiXa 148
## 3 discordant aPTT above antiXa 42
## 4 high aPTT low antiXa 1
## 5 low aPTT high antiXa 1## # A tibble: 5 × 2
## # Groups: concordance_low_simple [5]
## concordance_low_simple n
## <fct> <int>
## 1 concordant 327
## 2 discordant aPTT below antiXa 183
## 3 discordant aPTT above antiXa 59
## 4 high aPTT low antiXa 1
## 5 low aPTT high antiXa 1##
## Counts (Classic):## # A tibble: 2 × 2
## classic_discordant n
## <lgl> <int>
## 1 FALSE 49
## 2 TRUE 63##
## Counts (Modeled):## # A tibble: 2 × 2
## modeled_discordant n
## <lgl> <int>
## 1 FALSE 77
## 2 TRUE 35##
## Classic vs Modeled Discordance Table:## modeled
## classic FALSE TRUE
## FALSE 43 6
## TRUE 34 29## Warning: 'tidy.table' is deprecated.
## Use 'tibble::as_tibble()' instead.
## See help("Deprecated")## # A tibble: 4 × 3
## classic modeled n
## <chr> <chr> <int>
## 1 FALSE FALSE 43
## 2 TRUE FALSE 34
## 3 FALSE TRUE 6
## 4 TRUE TRUE 29##
## McNemar Test:##
## McNemar's Chi-squared test with continuity correction
##
## data: discord_table
## McNemar's chi-squared = 18.225, df = 1, p-value = 1.963e-05## # A tibble: 1 × 4
## statistic p.value parameter method
## <dbl> <dbl> <dbl> <chr>
## 1 18.2 0.0000196 1 McNemar's Chi-squared test with continuity corr…## === 3×3 Classification Table (Classic vs Polynomial Model) ===## modeled
## classic aPTT < antiXa aPTT > antiXa concordant
## aPTT < antiXa 14 4 30
## aPTT > antiXa 0 11 4
## concordant 4 2 43##
## === Chi-square Test ===##
## Pearson's Chi-squared test with simulated p-value (based on 5000
## replicates)
##
## data: discord_table_full
## X-squared = 55.029, df = NA, p-value = 2e-04##
## === Discordant-Only 2×2 Table (Low vs High) ===## modeled
## classic aPTT < antiXa aPTT > antiXa concordant
## aPTT < antiXa 14 4 30
## aPTT > antiXa 0 11 4
## concordant 4 2 0## Warning: 'tidy.table' is deprecated.
## Use 'tibble::as_tibble()' instead.
## See help("Deprecated")## # A tibble: 9 × 3
## classic modeled n
## <chr> <chr> <int>
## 1 aPTT < antiXa aPTT < antiXa 14
## 2 aPTT > antiXa aPTT < antiXa 0
## 3 concordant aPTT < antiXa 4
## 4 aPTT < antiXa aPTT > antiXa 4
## 5 aPTT > antiXa aPTT > antiXa 11
## 6 concordant aPTT > antiXa 2
## 7 aPTT < antiXa concordant 30
## 8 aPTT > antiXa concordant 4
## 9 concordant concordant 0##
## === Chi-square Test (Discordant Types Only) ===##
## Pearson's Chi-squared test with simulated p-value (based on 5000
## replicates)
##
## data: discord_types
## X-squared = 33.891, df = NA, p-value = 2e-04## # A tibble: 1 × 4
## statistic p.value parameter method
## <dbl> <dbl> <lgl> <chr>
## 1 33.9 0.000200 NA "Pearson's Chi-squared test with simulated p-val…6. Do counts of complications by discordance type differ when using classical (range-based) discordance categorization versus those identified with polynomial model?
Note that this analysis eliminated 1 lab results in which aPTT was above typical goal rangs and time-matched anti-Xa was below typical goal range and 1 lab result in which aPTT was below typical goal range and time-matched anti-Xa was above typical goal range.
##
## =============================
## Category: aPTT < antiXa
## =============================
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 48 20 0.417
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 18 4 0.222
##
## 2×2 Thrombosis Table:
## thrombosis no_thrombosis
## classic 20 28
## modeled 4 14
##
## Fisher Test:
##
## Fisher's Exact Test for Count Data
##
## data: table_mat
## p-value = 0.1649
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
## 0.6438088 11.8457795
## sample estimates:
## odds ratio
## 2.467387
##
##
##
##
## =============================
## Category: aPTT > antiXa
## =============================
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 15 6 0.4
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 17 6 0.353
##
## 2×2 Thrombosis Table:
## thrombosis no_thrombosis
## classic 6 9
## modeled 6 11
##
## Fisher Test:
##
## Fisher's Exact Test for Count Data
##
## data: table_mat
## p-value = 1
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
## 0.2311049 6.4469852
## sample estimates:
## odds ratio
## 1.214552
##
##
##
##
## =============================
## Category: concordant
## =============================
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 49 18 0.367
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 77 34 0.442
##
## 2×2 Thrombosis Table:
## thrombosis no_thrombosis
## classic 18 31
## modeled 34 43
##
## Fisher Test:
##
## Fisher's Exact Test for Count Data
##
## data: table_mat
## p-value = 0.4607
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
## 0.3280859 1.6249909
## sample estimates:
## odds ratio
## 0.7361483##
## =============================
## Category: aPTT < antiXa
## =============================
## # A tibble: 1 × 3
## n hemorrhage_n hemorrhage_rate
## <int> <int> <dbl>
## 1 48 2 0.0417
## # A tibble: 1 × 3
## n hemorrhage_n hemorrhage_rate
## <int> <int> <dbl>
## 1 18 1 0.0556
##
## 2×2 Hemorrhage Table:
## hemorrhage no_hemorrhage
## classic 2 46
## modeled 1 17
##
## Fisher Test:
##
## Fisher's Exact Test for Count Data
##
## data: table_mat
## p-value = 1
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
## 0.03641973 46.14255757
## sample estimates:
## odds ratio
## 0.74273
##
##
##
##
## =============================
## Category: aPTT > antiXa
## =============================
## # A tibble: 1 × 3
## n hemorrhage_n hemorrhage_rate
## <int> <int> <dbl>
## 1 15 3 0.2
## # A tibble: 1 × 3
## n hemorrhage_n hemorrhage_rate
## <int> <int> <dbl>
## 1 17 4 0.235
##
## 2×2 Hemorrhage Table:
## hemorrhage no_hemorrhage
## classic 3 12
## modeled 4 13
##
## Fisher Test:
##
## Fisher's Exact Test for Count Data
##
## data: table_mat
## p-value = 1
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
## 0.09859548 5.98490546
## sample estimates:
## odds ratio
## 0.817771
##
##
##
##
## =============================
## Category: concordant
## =============================
## # A tibble: 1 × 3
## n hemorrhage_n hemorrhage_rate
## <int> <int> <dbl>
## 1 49 3 0.0612
## # A tibble: 1 × 3
## n hemorrhage_n hemorrhage_rate
## <int> <int> <dbl>
## 1 77 3 0.0390
##
## 2×2 Hemorrhage Table:
## hemorrhage no_hemorrhage
## classic 3 46
## modeled 3 74
##
## Fisher Test:
##
## Fisher's Exact Test for Count Data
##
## data: table_mat
## p-value = 0.6769
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
## 0.2058079 12.4748253
## sample estimates:
## odds ratio
## 1.602364##
## =============================
## Low Range vs Polynomial — Thrombosis
## Category: aPTT < antiXa
## =============================
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 65 24 0.369
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 18 4 0.222
##
## 2×2 Table:
## thrombosis no_thrombosis
## low_range 24 41
## modeled 4 14
##
## Fisher Test:
##
## Fisher's Exact Test for Count Data
##
## data: table_mat
## p-value = 0.2766
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
## 0.5519251 9.4577870
## sample estimates:
## odds ratio
## 2.032359
##
##
##
##
## =============================
## Low Range vs Polynomial — Thrombosis
## Category: aPTT > antiXa
## =============================
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 19 8 0.421
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 17 6 0.353
##
## 2×2 Table:
## thrombosis no_thrombosis
## low_range 8 11
## modeled 6 11
##
## Fisher Test:
##
## Fisher's Exact Test for Count Data
##
## data: table_mat
## p-value = 0.7419
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
## 0.2855128 6.3885568
## sample estimates:
## odds ratio
## 1.322697
##
##
##
##
## =============================
## Low Range vs Polynomial — Thrombosis
## Category: concordant
## =============================
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 28 12 0.429
## # A tibble: 1 × 3
## n thromb_n thromb_rate
## <int> <int> <dbl>
## 1 77 34 0.442
##
## 2×2 Table:
## thrombosis no_thrombosis
## low_range 12 16
## modeled 34 43
##
## Fisher Test:
##
## Fisher's Exact Test for Count Data
##
## data: table_mat
## p-value = 1
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
## 0.3573489 2.4714823
## sample estimates:
## odds ratio
## 0.9489871CONCLUSIONS
thrombotic and hemorrhagic events are common among ECMO patients
thrombotic events occurred frequently, despite heparin use in this sample
there is moderate correlation between aPTT and antiXa among heparinized ECMO patients
categorizing discordance based on aPTT range:anti-Xa range may not effectively recognize values which are relatively unusual (to each other) at CUIMC. The counts of patients who are discordant based on classic groupings based on usual CUIMC systemic anticoagulation goals are different (numerically higher) than those with relative comparison via polynomial model used in this approach. Our approach may better recognize truly unusual values.
initiation of ECMO induces a hypercoagulable state which is associated with thrombosis development
aPTT < antiXa is associated with increased risk of thrombosis compared to concordant and/or high aPTT > antiXa states (all patients with this state developed thrombosis in our study).
Likewise, aPTT < antiXa may be associated with a decreased risk of hemorrhage (non-significant)
underpowered for survival analysis
the relative aPTT:antiXa value is more important than absolute values. Classifying discordance by range alone (i.e. within wide goal ranges) may not appropriately recognize hyper- or hypocoagulable states.
Comments on discrepancy between logistic and survival analyses
Logistic regression treats discordance as fixed (ever vs never) or at first discordance. Cox regression treats it as time-varying, so:
If a patient is discordant late, their earlier follow-up counts as “non-exposed” in Cox.
Logistic regression collapses all follow-up into a single binary exposure.
Even with post-discordance events only, this leads to different denominators and risk attribution.
We are likely very underpowered given infrequent events and small discordant subgroups (n = 11 and n = 17) and a time-varying exposure.